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Abstract 

The energy distribution and type of the particle with the highest momentum in quark jets 
are determined for each of the five quark flavours making only minimal model assumptions. 
The analysis is based on a large statistics sample of hadronic Z° decays collected with 
the OPAL detector at the LEP e"'"e~ collider. These results provide a basis for future 
studies of light flavour production at other centre-of-mass energies. We use our results to 
study the hadronisation mechanism in light flavour jets and compare the data to the QCD 
models JETSET and HERWIG. Within the JETSET model we also directly determine 
the suppression of strange quarks to be 

% = 0.422 ± 0.049(stat.) ±0.059(syst.) 

by comparing the production of charged and neutral kaons in strange and non-strange 
light quark events. Finally we study the features of baryon production. 
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1 Introduction 



By identifying the flavour of the quark from which a jet develops one can experimentally 
test both electroweak and QCD theories. The power of flavour tagging has been 
demonstrated in many studies of bottom and charm quark production. Tagging light 
quark jets is experimentally much more difficult as these jets are not as distinctive as 
bottom quark or charm quark jets. The main reason is that, unlike the heavy bottom and 
charm quarks, production of extra secondary up, down and strange quarks is abundant 
in jet development, making the identification of the hadron containing the primary quark 
ambiguous. Due to these difficulties, tagging of individual light quarks has been studied 
and used in only a few analyses, for example in . 

Whereas most of these analyses make assumptions about the details of hadronisation 
models, a method has been suggested in which reduces the reliance on these 
assumptions. This method has already been applied to determine the electroweak 
observables of individual light flavours by OPAL at the e^e~ collider LEP. In the 
present analysis, high-energy vr^, K^, Kg mesons, protons and A baryons are identified 
in the large Z" data sample and used as tagging particles. In addition, high-momentum 
e^, n^, D*^ mesons and identified bottom events are used to provide information about 
the heavy flavour backgrounds in these samples. As suggested in [^] and first confirmed 
by TASSO and more precisely studied in recent analyses, for example by SLD 0, 
these high-energy particles carry information about the original quark. In this paper we 
extend the method used in Reference to determine the probabilities i]l^{xp) for a quark 
flavour q to develop into a jet in which the particle with the largest scaled momentum 
Xp = 2pi/ a/s is of type i. 

The large number of Z° decays collected at LEP and their well-known properties give 
a unique opportunity for determining the probabilities ?7*(xp). From these measurements, 
we infer for the first time the flavour dependent fragmentation functions of light quarks. 
This allows us to study the hadronisation mechanism at an unprecedented level of detail. 
From these studies we determine in a direct way the suppression of strange quarks in the 
QCD sea and obtain insight into baryon production. 

Apart from such hadronisation studies, the results may also be applied to different 
environments. After taking into account QCD scaling violations which can be rather 
precisely determined, the rj^^Xp) from the Z° allow one to calculate the rj^^Xp) at other 
centre-of-mass energies. Possible applications include studies of light flavour production 
rates at other centre-of-mass energies |^ and the decay properties of the W boson, top 
quark or, if discovered, the Higgs boson. 

Section ^ contains a summary of the method. Section ^ describes the relevant features 
of the OPAL detector. The event selection and the tagging particle identification are 
described in Section ||. The determination of the rfq^Xp) is described in Section |] and 
their systematic uncertainties in Section ^. The results are shown and used to determine 
some properties of hadronisation in Section ^ 
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2 Method 



As detailed in the rf^i^Xp) are determined by using tags in event hemispheres^. Each 
event is separated into two hemispheres using the plane perpendicular to the thrust axis 
containing the interaction point. Each hemisphere is searched for the highest momentum 
particle, labelled i, subject to a minimum Xp requirement. If there are a number Nq of 
hemispheres which originate from a quark of type q and a number Nq^i{xp) of tagging 
particles i with a scaled momentum of at least Xp = Xcnt in these hemispheres, then the 
probability to find a tagging particle i with a scaled momentum of at least Xp is: 



Nq. 



The determination of the true viq{xp) at the "generator level", i.e. corrected for detector 
efficiencies and misassignment of the several particle types, is the main experimental 
aim of this paper. The particles considered are those which have a high probability to 
tag light flavours: vr^, K^, Kg mesons, protons, and A baryons. Charge conjugation is 
implied throughout this paper. What can be measured at the "detector level", i.e. before 
corrections for detector efficiencies etc., are the number of hemispheres tagged by a particle 
of type z, labelled Aj and called "single-tagged hemispheres", and the number of events 
containing a tagging particle in both hemispheres, labelled Nij and called "double-tagged 
events" , where i and j are the tagging particle types. 
These numbers are related to the probabilities: 

AT, 



AT, 



had 



i^v) = 2 Yl rtq'''^{xp)Rq (1) 

g=d,u,s,c,b 

and ^{xp) = (2-5,,) Y: P..(^p)<^^'(^p)^r''(^p)^'?' (2) 



had q=d,u,s,c,b 







where Sij = 1 if i = j and zero otherwise and A^had is the number of hadronic Z' 
decays. The superscript 'exp' denotes that the rj^^Xp) include possible distortions due 
to detector effects. The parameters pij{xp) take into account correlations between the 
tagging probabilities in opposite hemispheres, due to kinematic or geometrical effects, for 
example, and will not be equal to unity if such correlations exist. Rg is the hadronic 
branching fraction of the Z° to quarks q: 

- -F • 

J- had 

Rc and R^ are fixed to the LEP average measurements 0]. Given the good agreement of 
the Standard Model with data 0, in particular the agreement of the measured Rq, we 
fix -Rd/-Rught, -Ru/-Rught and -Rs/-Rught to their predicted values [0, such that J2qRq = I5 
where -Rught = Rd + Ru + Rs- 

^In this analysis, wc denote hemispheres as representing quark jets, since we are interested in studying 
the evolution of primary quarks into different hadron types. 
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The true f]g{xp) are found after correcting for detector efficiencies and misassignment of 
the tagged samples. The relationship between the true r]^g{xp) and the observed 77^' '^^'^{xp) 
is parametrised by a flow matrix, S^, which is taken from the simulation: 

^r^'(^p) = E^Hi^p) (3) 

i 

N ■ (x 

Ng_.,{xpr^ ' ^ ' 

where the sum over i includes all tagging particle types at the generator level and 
Nq^i^j[xp)^'" is the number of g-flavour Monte Carlo hemispheres tagged by particle i at 
the generator level but j in the detector. Sj is found to vary slowly with Xp. In addition 
it is necessary to count events which are untagged at the generator level but still give rise 
to a tagging particle in the detector and will henceforth be denoted "other background" . 
For example, these events can be tagged by a particle which is not considered in this 
analysis, such as a stable hyperon (e.g. or are tagged by a high- momentum 

particle which has a true momentum slightly below the minimum required Xp due to the 
finite momentum resolution of the detector. 

The system of equations (|^) and (||) has 20 equations (5 single and 15 double tags) 
with 25 unknown viq{xp) for the five quark fiavours produced in TP decays. We extend the 
system of equations in the following two ways: 

1. In order to better measure the heavy fiavour i]q{xp), we include charm and bottom 
tags by identifying D*^ mesons, or a vertex displaced from the interaction point. 
These techniques have been used previously in OPAL papers |jlT|,|T^ and are briefiy 
described in Section |4.4| . Charged leptons e^ and /i^, which mainly tag heavy 
fiavours but are still a source of background in the light fiavour charged hadron 
samples, are also identified and included in the equation system. Note that the 
vertex tag does not depend on Xp. 

2. In order to reduce the number of unknown r]q{xp), we invoke hadronisation 
symmetries such as 'qj'^ = r]'^^ , which are motivated by the fiavour independence of 
QCD and SU(2) isospin symmetries. They have been extensively discussed in 0] 
for Xp > 0.5. At lower momenta the relations are potentially broken by isospin 
violating decays, for example 0(1020) to charged and neutral kaons. Nevertheless, 
the relations 



Vd = Vn , 

vf^ = vf and 

e± e± 
= Vn 

are expected to be valid to high precision also after decays. This has been checked 



using the QCD model JETSET |13| after adjusting the yield and energy dependence 



of prominent resonances to the measurements at LEP |14|. We find the relations 



hold to within 2% above Xp = 0.2, the range used in this analysis. Here K° is made 
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up of both Kg and K^, which are assumed to be equal. A relation that is used which 
is violated by up to 10% at low Xp ~ 0.2 due to decays is 

„A(A) _ „A(A) 

When introducing these hadronisation symmetries into the equation system, we 
make whatever small corrections for isospin-violating decays are necessary according 
to the JETSET Monte Carlo. The HERWIG model H is not used to check the 
hadronisation symmetries because it violates SU(2) isospin symmetry for technical 
reasons [0. 

These additions give a total of 54 equations with 41 unknown r]^^{xp). The equations 
are solved requiring a minimum Xp > 0.2, 0.3, 0.4, 0.5 and 0.6. In this paper, full results 
are presented for a minimum Xp > 0.2 and summarised for the other cut values. 



3 The OPAL Detector 



The OPAL detector is described in detail in The relevant features for this analysis 
are summarised in this Section. OPAL uses a right-handed coordinate system, where the 
2;-axis points along the electron beam, r is the coordinate normal to this axis, and 6 and 
are the polar and azimuthal angles with respect to z. 

The central tracking system, inside a 0.435 T axial magnetic field, provides a charged 
track momentum resolution of CTp/p = 0.02 © 0.0015 pt, where pt is the momentum 
component perpendicular to the beam axis in GeV. A silicon microvertex detector ||18|| , 
close to the interaction point, is surrounded by three drift chambers: a vertex detector, 
a large volume jet chamber which provides up to 159 space points per track, and z- 
chambers which give a precise measurement of the polar angle of charged tracks. The 
large number of samplings in the jet chamber also provides a determination of the specific 



ionisation energy loss, dE/dx, with a resolution of a{dE/dx)/{dE/dx) ~ 0.032 |]19[ in 
multihadronic events for tracks with | cos^| < 0.7 and the maximum number of samplings. 
At larger | cos 6*1 the resolution is degraded because fewer measured points are available. 
The dE/dx measurements have been calibrated using almost pure control samples of, 
for example, pions from Kg, //-pair events and photon conversions into electrons, such 
that the central values are known to 0.10a{dE/dx) and the resolution to a precision 
of 10% [EDI. The electromagnetic calorimeter consists of 11704 lead glass blocks, each 



subtending a solid angle of 40 x 40mrad^. The muon chambers surround the calorimeter, 
behind approximately eight absorption lengths of material. 

Detector efficiencies and possible detector biases are studied with approximately six 
million simulated hadronic decays generated with the JETSET 7.4 model [T^ and 
passed through a detailed simulation of the OPAL detector ^1||. The fragmentation 
parameters have been tuned to describe event shapes and other distributions as described 
22[. In addition, for fragmentation studies one million fully-simulated hadronic Z° 



m 



decays generated with the HERWIG 5.8 Monte Carlo generator are used. 
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4 Event Selection and Tagging Methods 



The analysis uses approximately 4.1 million multiliadronic T? decays collected between 
1991 and 1995. The standard OPAL multihadronic selection is applied To select 

events which are well contained in the detector, the polar angle of the thrust axis, ^t, 
calculated using charged tracks and electromagnetic calorimeter clusters which have no 
associated track in the jet chamber is required to satisfy | cos6't| < 0.8. To assure good 
bottom quark tagging quality it was also required that the silicon microvertex detector 
be functioning well. 

In this analysis we select tagging particles with Xp > 0.2. The selection of these 
highly energetic particles enhances the background fraction from t^t~ events, so 



we require in addition each event to have at least eight well- measured tracks [24|. With 



these requirements, 2 820 220 events are retained. In this event sample the r background 
is reduced to less than 0.03 %, as estimated using fully simulated events generated with 
the KORALZ Monte Carlo generator and so can be neglected. 

Next a high-energy stable particle or a charm or bottom tag is required. The selection 
of particles was optimised for the highest accuracy of the desired 77* (xp), balancing the 
potential loss in separation power against efficiencies. As discussed in Section |], we look 
for the particle (vr^, K±, p(p). Kg, A (A), e±, /i^ or D*^, ) in each event hemisphere 
with the highest scaled momentum Xp. To ensure good charged pion, kaon and proton 
separation, and reliable Kg and A reconstruction, we require that the tagging particles 
have polar angles | cos6'| < 0.9. 



4.1 Stable hadrons 

The dE/dx measurement of good quality tracks is used to identify charged pions, charged 
kaons and protons. For each track the dE/dx weight Wi is used to separate the particle 
types. The weight is the probability that the track is consistent with a hypothesised 
particle i. We require: 

• for pion candidates: > 0.1 and Wk± < 0.1; 

• for kaon candidates: Wk± > 0.1 and Wj^i < 0.1; 

• for proton candidates: Wp(p) > 0.1 and Wk± < 0.1. 

These selection criteria give three disjoint samples. 
Averaged over all five quark fiavours, 

express the probability that a particle identified as type j stems from a true particle type 
i. The values are given in Table |1| for samples with Xp > 0.2. The determination of 
the fiow matrix Sj, which is taken from simulation, was discussed in Section |^. As also 
discussed in Section 0, a few percent of the tagging particles have a true Xp value below 
the cut imposed on the measured Xp value but are tagged in the detector due to the finite 
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momentum resolution. The Monte Carlo also predicts that there is a 1% contamination 
from charged hyperons, mostly S~, in the proton sample. These sources of background 
are included in "other background" given in Table |I|. The cuts, including the event and 
thrust axis cut, lead to the efficiencies shown in the bottom row of Table 0, which are 
defined as the number of hemispheres which are correctly tagged at the detector level 
divided by the number of hemispheres which are tagged at the generator level. 



4.2 Electron and muon identification 

Electrons are identified using a number of discriminating variables, principally the dE/dx 

Muons are selected by 
As can 
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and the energy loss in the electromagnetic calorimeter 
matching tracks in the central detector with hits in the muon chambers 
be seen from Table |l|, the efficiencies to correctly tag an electron or a muon are about 
20% and 70% in this hadronic jet environment, respectively, with purities of around 60%. 



4.3 Kg and A identification 

The procedures to identify the weakly decaying Kg and A are described in |^ and |2l 
respectively. The decays Kg — > tt+tt" and A — > pvr" are reconstructed by combining 
two oppositely-charged tracks which have a crossing point in the plane orthogonal to 
the beam axis. If a secondary vertex is found, the invariant masses and mp^- of 

the 7r+7r~ and pvr" mass assignments are calculated. Kg candidates are required to have 
invariant masses in the ranges 430 MeV < < 570 MeV and mp^- > 1.13 GeV, to 

reduce the contamination from A pvr" decays. Similarly, all candidates which have 
1.10 GeV < mpTT- < 1.13 GeV are accepted as A candidates. The Kg selection in ||2^ is 
extended in the present analysis to | cos6\ < 0.9 from 0.7, resulting in a slightly worse 
overall mass resolution, but the acceptance is increased and is the same as the other 
particle tags used, thus reducing geometric hemisphere correlations. The combinatoric 
backgrounds are estimated from the Monte Carlo and cross-checked by determining the 
backgrounds from candidates with invariant masses in sidebands around the signal. 



4.4 Charm quark and bottom quark tags 

The sample enriched in charm quark events is found by selecting hemispheres with a 
high-energy D*^ [|l^. The decay modes and the cuts used on the x^* values, which are 
calculated from the measured decay products of the D*"^ candidate, are 

D*+ 



D07r+ 










U K- 


'n~ 


h 




> 0.4 


U K- 


-e^ 




> 0.4 


U K- 


>" 






> 0.4 


U K- 








> 0.4 


U K- 








> 0.5 



In the simulation the D*^ tagging efficiency is about 1% for a charm quark purity of about 
58%. 
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The bottom quark tag uses a number of discriminating variables calculated from a 
reconstructed secondary vertex [0 . From the number of double and single-tagged events 
the hemisphere tagging efficiency is found to be about 19% for a bottom jet purity of 
about 96%. These efficiencies and purities are used only for the cross-check outlined in 
Section |6]^. The hemispheres tagged as bottom are counted even if they are already 
present in the high Xp tagged samples. 



5 Determination of Tj^q{xp) 



The numbers Ni and Nij of measured single- and double-tagged events are given in Table ^ 
for Xp > 0.2. They are used as input to the equation system which is solved for the 77* by 
using a fit- The function is defined as: 



(6) 



where 



Ni = N,-J2{1 + 5,,)N,, and 
j 



j 



are used to correct for double- counting of hemispheres in the single- and double-tagged 
samples. These two corrections are necessary to remove double-tagged events from the 
sample of single-tagged hemispheres. 

In addition, the hadronisation symmetries given in Section ^ are used after being 
corrected for detector effects and making small corrections for isospin-violating decays 
according to the JETSET Monte Carlo. Furthermore, certain very small, and therefore 
unmeasurable, rj^ are fixed to their JETSET values at the generator level, namely: ?7d ^ .,5 
rj^l^s and i]^'^^'^- The pij parameters, parametrising possible kinematic and geometrical 
correlations, are taken from the simulation. Geometrical correlations lead in general to a 
positive correlation pij > 1. Motivated by simulation studies, the correlation is assumed 
to be the same for all tagging particle types except for the D*^ and the bottom tag. 
Typical values are pij = 1.020 ± 0.002 at Xp > 0.2 and = 1.13 ± 0.03 at Xp > 0.5, 
where the errors are from Monte Carlo statistics, and i and j run over all tagging particle 
types except for D*^ and the bottom vertex tag. For the D*^ and bottom vertex tags, 
the correlations are determined individually for each measured double-tagged sample. For 
example, p7r±D*± = 1-048 ± 0.021 and P7r±b-vtx = 1-018 ± 0.006 for Xp > 0.2, where the 
errors are again from Monte Carlo statistics. The extracted 77*' '^^^{xp) are corrected for 
the detector effects using the flow matrix Sj in equation 0. 

The results after corrections for detector efficiency and particle misassignment are 
listed in Table |^ for Xp > 0.2. The table also includes the statistical and systematic 
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uncertainties and a comparison with the JETSET and HERWIG models. We give details 
of the results only for the tagging particle types which mainly tag light flavours, namely 
vr^, K^, Kg, proton and A. The statistical correlations between the parameters are given 
in Table |. 

The corrected results also for Xp cuts other than Xp > 0.2 are summarised in Table ^ 
with statistical and systematic error combined. Some of the larger rj^gixp) are shown in 
Figures |I]- Correlations between the rjg{xp) for different particle types and between the 
values obtained with different Xp cuts are discussed in Section 0. 

The solutions were checked to be unique and that the error matrices were positive 
definite. The per degree of freedom of the solutions are typically ~ 1.2, and are given 
in Table |[ 



6 Systematic Uncertainties 

The validity of the method used in this paper is tested using approximately six million 
hadronic Z° decays generated using the JETSET Monte Carlo and including a full 
simulation of the OPAL detector. The rj^^Xp) obtained from solving the equation system 
agree with the Monte Carlo predictions. 



6.1 Main uncertainties 

The main sources of systematic uncertainty are due to the limited knowledge of the 
efficiencies and purities of the particle identification. Others are due to the flavour 
composition of the D*^ and bottom-tagged samples. The third class of uncertainties 
is related to opposite-hemisphere correlations in the double-tagged samples. Since these 
classes of errors are largely uncorrelated, we estimate their individual impact on the 
rj^q{xp) and add them quadratically to obtain the overall systematic error. The errors 
are determined by changing in turn each input parameter according to the estimated 
individual range of uncertainty, repeating the analysis, and interpreting the shifts as 
the error contribution. A break down of the individual error contributions for the most 
important rf^i^Xp) is listed in Table | for Xp > 0.2. Relative contributions to the systematic 
error at other minimum values of similar. 

The following systematic uncertainties are considered: 

• Charged particle purity and efficiency: 

Systematic errors are applied to the charged pion, charged kaon, and proton yields. 
The uncertainties in these corrections are estimated by varying the widths and mean 
values of the ionisation energy loss in the simulation according to the uncertainties 
discussed in Section |^ [^. These errors are the dominant ones for all rf ixp) of 



charged hadrons. 



The uncertainties in the electron and muon identiflcation have been discussed in |T2 . 
The error of the electron identiflcation efficiency is due primarily to uncertainties in 
modelling the dE/dx. The modelling of the muon efficiency has been checked using 
/i-pair and e"'"e~ — >• e"'"e~/i"'"/i~ events. The effects on the hadron f]g{xp) are small 
and are included in the error due to charged particle efficiency and purity. 



11 



Efficiencies of Kg and A: 

The uncertainties of tlie Kg and A efficiencies as given for Xp > 0.2 in Table 
for example, are taken into account. Since the relative yields of Kg and K^ are 
important for the separation of up and down quark jets, the uncertainty contributes 
significantly also to the 77^ , for example, to rj^ and 77^ . The relevant sources of 
systematic error are described in and . For the Kg the systematic errors for 



the region 0.7 < | cos^^| < 0.9 were taken to be double those in the barrel region, 
motivated by the factor of two worse mass resolution in the endcap. 

Charm tag efficiency: 

The relative uncertainty in the D*^ reconstruction efficiency was conservatively 
estimated to be ±10%. This source of error has a negligible effect on the results. 

Hemisphere correlations: 

Correlations due to kinematic and geometrical effects are accounted for by the pij 
parameters, which are taken from Monte Carlo simulation. The values of pij are most 
sensitive to changes in the angular acceptance of the tagging particles and the thrust 
angle cut. Variations in maximum | cos6'| of the tagging particles between 0.7-0.9, 
and different cuts on the maximum | cos^tI between 0.7-0.9 show that the changes 
of the Pij are well simulated. A ±0.01 absolute systematic error, representing the 
maximal disagreement between data and Monte Carlo, is assigned for the simulation 
of the Pij values. 

Other background: 

Contributions to the detector level r^* from events which are not tagged at the 
generator level are taken from the JETSET Monte Carlo events. Such events are 
mainly due to tags which have a true Xp lower than the minimum Xp cut used but are 
tagged due to the finite momentum resolution in the detector, spurious tracks, and 
combinatoric background in the case of the Kg and A samples. Another source of 
other background (mainly in the proton sample) is due to stable charged hyperons, 
mostly S". 

These backgrounds represent either an absolute contribution to 77*' '^^^{xp) or a 
constant background fraction which scales with the rig{xp). The systematic errors 
on the estimations of these backgrounds are taken as the differences in the viq{xp) 
if the analysis is repeated under the two assumptions, namely treating the other 
background events as a fraction of the detector level rf^ '^^^{xp) or as an absolute 
contribution. This procedure takes into account uncertainties in the JETSET 
modelling of both the magnitude and the Xp dependence of the background sources. 

Charm tag background: 

The flavour composition of the D*^ sample has been discussed in [O]. The 



fraction of bottom quark jets in the D*^ sample can be directly determined. The 
contribution from gluon splitting — >^ cc is negligible. The flavour composition of the 
combinatorial background is taken from Monte Carlo. We varied the contributions 
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individually by ±50%. The corresponding uncertainties have been included in the 
systematic errors. 

• Fixed quantities: 

The quantities which were fixed in the fit, namely us) VdTs^ i]]^'^''^^, are each 
in turn varied by ±100% and the corresponding shifts in the 77* taken as systematic 
errors. 

• Hadronisation symmetries: 

As discussed in Section ^ the hadronisation symmetries used to solve the equation 
system may be broken by up to 2% at low Xp and 10% for the relation between 
Vdi^p) Vui^p) 0)0- The relations are corrected for any breaking which is 
present in the JETSET Monte Carlo. Assuming a systematic error equal to the 
maximal allowed breaking, the rj^gixp) change only marginally. 

• Z° branching ratios: 

The uncertainties due to the branching ratios Rq into quarks have been estimated 
by varying each fraction within certain limits. In the case of bottom and charm 
quarks these are given by the rather precise measurements at LEP |p. Branching 
fractions into the individual light quarks are less well determined. A direct 
measurement has only been performed in 0, which we do not consider because it 
used a variation of the method applied in this paper. However, there are constraints 
on the electroweak couplings of up and down quarks coming from lepton nucleon 
scattering and final-state photon radiation from quarks which agree with the 
Standard Model expectation. 

To take into account uncertainties in the light fiavour Rq, we vary Ri/{Rd + Ru + Rs) 
where i = d, u, s by ±10% from the Standard Model values taking into account their 
well measured sum 1 — Ry, — Rc = 0.606 ± 0.010 g|. The rf^^Xp) values change by 
a maximum of 0.5%, and the only marginally. Since we assume the Standard 
Model in this analysis, we do not include this small source of error in the overall 
systematic errors. 

6.2 Cross-check on events with a heavy quark tag 

We followed the procedure detailed in in order to make a cross-check of the principal 
results. Compared to the light fiavour tags based on high Xp stable hadrons, the purity 
of heavy quark tags is much higher. The cross-check makes use of the charm and bottom 
tag efficiencies and purities from Monte Carlo, mentioned in Section [4.4| . By counting the 
number of light fiavour tags in event hemispheres opposite to a heavy fiavour tag, one can 
determine the r/*' '^^^{xp) directly without using a large system of equations. This method 
leads to results which are consistent with those from the main method used in this paper. 
The agreement between data and the JETSET model is found to be in general quite good. 
The biggest discrepancy between the data and the prediction of the JETSET model is for 
?7^(xp), which will be discussed in more detail in Section |7[1|. 
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7 Results and Hadronisation Studies 



The viq{xp) values with their statistical and systematic uncertainties are listed in Table |] 
for Xp > 0.2. The largest in light flavours for mesons are shown in Figures |1| and as 
a function of the cut on Xp. In addition, for baryons, the r]^, r^^, t]^ and r]^ are shown in 
Figures | and ^. 

In most cases only weak correlations exist between the r/* for different tagging particles 
i, though stronger correlations exist between different flavour 77* with the same tagging 
hadron. The statistical correlation coefficients for the most important r/* are given in 
Table ^ for Xp > 0.2. These correlations are typical also for the other minimum values of 

Xp. 

In all cases the expected pattern of the leading particles holds: the up and down 
quarks fragment mostly into pions, whereas the strange quarks fragment mostly into 
and Kg, although the fraction of vr^ is sizable also for strange quarks, especially at low 
Xp. The heavy charm and bottom quarks produce mostly high-energy pions and kaons. 



7.1 Comparison to JETSET and HERWIG 

We compare the dominant fragmentation functions for the different flavours to the 
expectations of the HERWIG and JETSET models. The results are shown in Figures ^ 
and ^ for mesons, and in Figures |^ and § for baryons. Note that the data points are 
shown for different minimum values of the Xp and therefore are correlated. For x > 0.2 
the JETSET and HERWIG expectations for all determined fragmentation functions are 
listed together with the data in Table ^ 

Hadronisation is quite differently modelled in the two QCD generators. Whereas 
JETSET uses the Lund string model |2^, HERWIG invokes principally the cluster decay 



mechanism [29|. Both models, JETSET more than HERWIG, contain several parameters 
which cannot be derived from first principles. For this comparison we use the standard 
OPAL tuning p2[ which is optimised to describe the overall event properties and inclusive 



particle production. Our measurements of the flavour dependence of the fragmentation 
function allows us to test the correctness of the model at a new level of detail. 

Most of the tagging probabilities are well reproduced by both JETSET and HERWIG. 
Exceptions are the consistent underestimation in HERWIG of the meson production in 
bottom events. In addition HERWIG seems to underestimate rj'^ and both HERWIG and 
JETSET seem to underestimate 77^ . The significance of these deviations is, however, 
only at the level of two standard deviations. The distributions for kaons have a similar 
shape for the various quark species. The distributions for pions are significantly steeper 
than those for kaons, which can at least partly be explained by the larger fraction of pions 
from resonance decays that will be found at lower Xp values. Particularly rj^, which can 
only be due to either decays or if both an up and a down quark are produced from the 
hadronisation sea, is steeper than rj^ and 

Whereas proton production is reasonably described by JETSET, the HERWIG 
prediction deviates significantly from the data. The flavour integrated rate pO| 



IS 



overestimated at high Xp in this model. As can be seen from Table |^, especially the 
fractions of protons in up and down quark events are higher. In fact the excess in the 
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HERWIG prediction is almost exclusively due to these quarks, which are in HERWIG for 
Xp > 0.2 about twice as high as in the data. The data also show that the yield of protons 
from up quarks is higher than that from down quarks. This will be discussed in more 
detail in Section |7^. 

HERWIG also significantly overestimates A production for all light quark species, as 
can be seen in Figure ^ The overestimation for ri^{x > 0.2) is less pronounced; however, 
the shape of the fragmentation function is softer than in the data. In the case of A baryons 
in strange and charm events, the JETSET expectation differs from the data as seen in 
Figure ^. The s — >^ A yield in the data is only about half of that expected although the 
shape is consistent. For c — A the yield is underestimated by a factor of 2—3 and the x- 
dependence tends to be steeper. To study whether the discrepancy in the rate may be due 
to the analysis procedure, we compare A production directly in data and in the JETSET 
simulation including detector effects. To enrich charm and strange events, respectively, 
we search for A in hemispheres opposite to a tagged D*~ or K+, and A production in 
hemispheres opposite to a tagged D*"^ or K~. The resulting pvr" mass spectra are shown 
in Figure ^ The underestimation of the A production in charm events in the simulation 
is clearly visible, as is the overestimation at high Xp of A production in strange events. 

In addition to studying absolute rates of individual particle species for a specific 
flavour, as the next step we compare relative yields for the same flavour or the same 
particle type. These relations may reveal symmetries in the hadronisation mechanism. 

7.2 Strange quark fraction in QCD vacuum 

In a next step we compare the yield of in up and strange quark events and Kg in down 
and strange events. Within JETSET the ratio of the production yields of the primary 
hadrons is a direct measure of 7s = V{s)/V{u,d), i.e. the relative quark production 
probabilities in the hadronisation sea. We present the results in Figure ^ and Table |^. 
The full lines show the expected ratio in JETSET for a 7^ value indicated by the dotted 
lines. The difference of up to 10% between the expected ratio and 7s is due to decays, 
particularly of the L = 1 meson supermultiplet. The comparisons show that i]^^ /ri^"^ 
(Figure |^a) and rj^ /rjf (Figure ^d) are good estimators of 7s. 

No significant dependence on Xp is observed for either the or Kg measurements, 
which is consistent with expectations. A combined K^ and Kg analysis is made by 
invoking SU(2) isospin symmetry which implies that r]^"^ = 1]^° . After taking into account 
correlations and correcting for isospin-violating decays at Xp > 0.2, we obtain 

7s = 0.422 ± 0.049(stat.) ±0.059(syst.) 

The systematic uncertainty includes an error of 0.042 to take into account variations of 
the correction factors due to the uncertain amounts of resonance production, found by 
varying the contributions of the L = 1 meson supermultiplet by ±50%. 

This value of 7s is consistent with, although somewhat larger than, previous 
measurements [^] which are, however, in most cases rather indirect. Comparing the 
data in more detail with the JETSET prediction in Table ^ one observes that the rj^ 

and ?7s ^ are in good agreement. However, in the data more K are found in up quark 
events and more Kg in down quark events than predicted by JETSET. 
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In the case of HERWIG the ratios of the flavour dependent Kg and K production 
have no simple interpretation in terms of a single parameter. The ratios are fairly similar 
to those of JETSET and thus agree with the data. 



7.3 Baryon hadronisation 

The mechanism of how three quarks coalesce in the jet development to form a baryon is 
still a puzzle of hadronisation. Our measurement of the flavour dependence of the proton 
and A yields provides additional new input. 

The ratio rj^/rj^ is shownQ in Figure ^ and listed in Table ^. Within the LUND string 
model ideally the ratio Tj^/rj^ at high Xp would be a direct measure of the size of the 
suppression of diquarks |]32| with spin 1 relative to spin 0, since Fermi statistics requires 
a (uu) diquark to have angular momentum L = 1. However, decays from heavier baryons 
such as A or A resonances tend to change the ratio. Although our result agrees with the 



production of diquarks as suggested in [32| and already supported by studies of baryon 



number compensation in jets the uncertainties are so large that the data are also 
consistent with models that form baryons from quarks that are statistically produced in 
rapidity. HERWIG, which incorporates a democratic production of diquarks, predicts the 
production of protons from u and d jets to be more equal than JETSET. The data tend 
to be smaller than the HERWIG expectation. 

We also observe the suppression of strangeness in baryon production by measuring the 
ratio of A baryon production in down and strange quark events. Note that in solving the 
equation system we have assumed that rj^ ~ rj^. After the production of the primary 
down or strange quark, A baryons are formed by picking up a pair of (us) or (ud) quarks, 
respectively. Indeed fewer A baryons are found in down (and hence up) quark jets than 
in strange jets. The suppression agrees with the JETSET and HERWIG models. The 
results are shown in Figure |^ and listed in Table ^. 

Finally we compare the production of baryons and mesons in events of the same 
primary quark type. The ratio of proton to pion production in up quark events and 
the ratio of A to charged kaon production in strange events are given in Table ^ and 
shown in Figures ||a and §b, respectively. The JETSET expectations fall above the 
measured data points. Although these measured ratios are poor estimators of the level 
of diquark suppression (indicated by the dotted line) within the JETSET model due 
to large contributions from decays, the suppression of baryons relative to mesons is 
clearly observed. For both ratios the HERWIG expectation is significantly above the 
measurement as already mentioned in Section |7.1| . 

In addition one can form the double ratio, {vt /vf)/{Vn/Vu) which within the JETSET 
model should measure the same quantity P(ud)/P(u), modified only by decays, where 
V{x) indicates the probability to pick out either a quark or diquark x from the QCD 
vacuum. The data yield for Xp > 0.2 



1.23 ±0.31, 



^The results for Xp > 0.30 are not very precise due to a lack of separation power between up and down 
quarks. 



16 



consistent with the JETSET expectation of 1.55, but significantly lower than the 
HERWIG prediction of 2.24. This indicates that the inclusive baryon production is 
badly modelled and also the relations between different meson or baryon species are 
unsatisfactorily simulated in HERWIG. 

8 Conclusions 

In this paper we have reported on a determination of the probabilities rfq{xp) of leading 
particles to originate from individual quark flavours in TP decays. We studied the 
production of leading n"^, K^, Kg, proton, and A for Xp > 0.2 up to 0.5. The measurement 
has only a minimal reliance on hadronisation models. In general we observe the expected 
behaviour that the flavour of the primary quark is reflected in the leading particle, i.e. up 
and down quarks lead mainly to highly energetic pions, while strange quarks lead mainly 
to kaons. 

These measurements allow several aspects of hadronisation to be studied rather 
directly, in contrast to many previous analyses which rely strongly on a model unfolding of 
different contributions. In particular we determine from the relative production of leading 
charged kaons in up and strange quark jets and leading Kg in down and strange quark 
jets the suppression of strange quarks in the QCD vacuum: 

73 = 0.422 d= 0.049(stat.) ±0.059(syst.) 

We also find that leading protons are more frequent in up than in down quark jets. We 
also observe the suppression of strange diquarks in (d, u) ^ A events and baryons relative 
to mesons in events of the same quark flavour. 

For most quark flavours and particle types the JETSET model reproduces the 
measurements well. A possible exception is the production of A baryons in charm quark 
events which appears to have a higher yield and a harder fragmentation function than 
expected. HERWIG provides in general a good description of mesons in light quark events 
but has deflciencies in baryon production, in particular the relative yields of different 
baryon types and the ratios of baryons and mesons in the same flavour jet. 

In addition to these hadronisation studies, our measurements of the r]g{xp) may also be 
interesting for future experiments. In providing tagging probabilities for light flavours with 
hardly any reliance on hadronisation models, the Ti^g{xp) can be applied at other centre- 
of-mass energies or in the study of heavy particle decays. This allows a determination 
of the light flavour production yields and properties in a model-independent way also for 
environments other than the Z°. 
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Assigned 


7r± 


K± 


P(p) 


c± 


True 




A(A) 


j3*± 


otlier 




0.790 


0.062 


0.003 


0.013 


0.007 


0.038 


0.005 


0.062 


0.019 


K± 


0.146 


0.568 


0.148 


0.002 


0.002 


0.017 


0.026 


0.071 


0.020 


P(P) 


0.040 


0.246 


0.551 


0.002 


0.001 


0.014 


0.081 


0.036 


0.031 


e± 


0.186 


0.023 


0.002 


0.620 


0.000 


0.024 


0.006 


0.128 


0.011 




0.100 


0.061 


0.002 


0.002 


0.643 


0.017 


0.007 


0.153 


0.015 




0.081 


0.030 


0.007 


0.004 


0.001 


0.691 


0.026 


0.101 


0.060 


A(A) 


0.047 


0.024 


0.024 


0.003 


0.001 


0.128 


0.696 


0.032 


0.045 


D*± 


0.143 


0.074 


0.019 


0.007 


0.006 


0.016 


0.012 


0.699 


0.024 


efficiency 


0.487 


0.441 


0.292 


0.228 


0.702 


0.155 


0.135 


0.033 





Table 1: Fractional compositions of the identified samples (rows) in terms of the true 
tagging particle, for Xp > 0.2. The dominant component of other tagged events (last 
column) is tagging particles which pass the minimum Xp requirement in the detector but 
whose true momenta are lower. The sum of the elements in each row is one. The last 
row gives the average efficiency to correctly tag a hemisphere, as taken from Monte Carlo 
simulation. Errors are discussed in the text. 
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Particle 
type 


Tagged 

1 • 1_ 

nemispneres 






I 

P(P) 


Double- 

e=^ 


tagged 


events 

T^O 

Ks 


A(A) 


D 


b-vtx 


7r± 


855043 
850442 


71601 
72515 


77576 
71411 


16287 
16258 


7127 
7146 


7689 
6637 


10425 
9299 


4717 
4404 


2483 
2238 


25074 
23221 


K± 


506538 
474123 




25717 
22 (42 


10120 

moo 

9182 


4135 

one o 

3988 


4475 

/I no 1 

4081 


7248 
6497 


3193 
3119 


1578 
1526 


14376 
13344 


P(P) 


101415 

100046 






963 
1019 


789 
839 


887 
815 


1375 
1219 


591 

583 


314 
290 


2744 

2784 


e± 


54370 
56235 








501 
479 


1294 
1186 


594 
542 


235 
253 


225 
219 


6219 
6221 


/^^ 


65029 
60767 










905 
838 


674 
541 


278 
274 


293 
262 


8898 
8254 


Kg 


71218 
64290 












523 
423 


454 
440 


239 

194 


2074 
1826 


A(A) 


31721 
30676 














107 
95 


111 
82 


1026 
973 


D*± 


17432 
16692 
















57 
76 


805 
791 


b-vtx 


245451 
246766 


















22472 
23047 



Table 2: Number of tagged event hemispheres and double-tagged events for Xp > 0.2. 
The upper numbers are for data and the lower for Monte Carlo, normalised to the same 
numbers of Z ^ gg events. 
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DPAT, Hafa 

WJTiilJ LIojLoj 


TF.TSF.T 

J 11/ ± oil/ J- 


llll/lXVV Ivjr 




U.c)oOD ± U.UUzr ± U.Uzo/ 


n ono^^ 
O.d9zo 


n 
O.dooo 




n 0001 _L n nnoi^ _l n noc;^^ 
U.oooi ± U.UUzD ± U.Uzoo 


n ooni 
0.0091 


0.0000 




U.i/Ul ± U.UUDZ ± U.Uioy 


n 1 OO y( 

0.1oo4 


O.io6/ 


Vc 


n 1 798 -1- n nnQ7 -i- n ni 


n 1 f^ns 

U.lOUo 


n 1 Asn 

U.140U 


Vh 




U. X^ZiU 


1 19Q 


Vd 


n ^ _i_ n ni no _i_ n nnon 
U.Uoi/ ± U.UiUz ± U.UUoU 


n n 1 ^ 
O.Ool / 


n n /I 1 
0.04ol 




[j.iZZi ± U.UioO ± U.Uz44 


O.O60 ( 


n n^no 
0.0 (06 




U.zoyU ± U.UUOO ± U.Uio/ 


n oon /I 


n oonn 
U.zoUU 


Vc 


u.uyoz zt U.UiUU zt U.UZUo 


nil 9^? 

U. i iZO 


nil n7 

U. i iU ( 


Vh 


U.UOZO ZE U.UUi( ZlZ U.UUOO 


u.uoou 


U.U4i04i 


p 
% 


0.0362 zt 0.0063 zt 0.0075 


n no 

0.0356 


n 0*70 r 

0.0795 




0.0569 zt 0.0086 zt 0.0109 


0.0666 


0.0913 




0.0328 ± 0.0042 ± 0.0084 


0.0232 


0.0319 


Vc 


n no/i^ _i_ n nnc^Q _l n nn^ri 
U.Uz4o ± U.UUoo ± U.UU/i 


U.UzOO 


n noQ /I 


p 


0.0232 ± 0.0010 ± 0.0045 


0.0253 


0.0199 


'Id 




0350 


0345 


l/U 


0228 zt 0107 zt 0100 


0.0251 


0229 


Kg 

71^ 


1210 zt 0028 zt 0096 


1161 


1160 


Kg 


0.0586 zt 0.0072 zt 0.0132 


0.0457 


0.0505 


K^ 

Vh 


n 097'? -1- n nni 9 -i- n nn99 

U.UZiO ZC U.UUiZ ZC u.uuzz 


n n99fi 
u.uzzo 


n ni Qc; 
u.uiyo 


A 
^d 


0.0231 ± 0.0025 ± 0.0020 


n n 1 ^0 

0.0172 


0.0566 


A 
Vn 


0.0211 ± 0.0023 ± 0.0020 


0.0158 


0.0542 


A 


0.0493 ± 0.0046 ± 0.0041 


0.0607 


0.1325 


Vc 


n n9Q'i -1- n nn7'i -i- n nnfi'? 

u.uzyo zn u.uufO zn u.uuoo 


n n9'ii 

U.UZO 1 


u.uoou 


'it 


0.0180 ± 0.0009 ± 0.001 1 


0.0182 


0.0170 



Table 3: Results for Xp > 0.2 after corrections for detector efficiency and particle 
misassignment. The first error shown is statistical and the second systematic. The two 
rightmost columns show the JETSET and HERWIG expectations with the OPAL tunings. 





vf 


vf 






vT 


p 
Vd 


Vl 


Kii 

Vd 






vf 


1.000 


-0.114 


-0.117 


-0.050 


-0.019 


0.023 


0.006 


-0.058 


-0.015 


-0.008 


vf 




1.000 


-0.013 


-0.118 


-0.334 


0.022 


-0.088 


0.081 


-0.029 


0.091 


vf- 






1.000 


-0.813 


0.250 


-0.618 


0.503 


-0.077 


-0.081 


0.067 


vf 








1.000 


-0.354 


0.541 


-0.603 


-0.009 


0.127 


-0.099 


vT 










1.000 


-0.091 


0.218 


0.168 


-0.206 


0.117 


p 

Vd 












1.000 


-0.825 


0.078 


-0.015 


0.086 


vl 














1.000 


-0.086 


0.047 


-0.013 


Vd 
















1.000 


-0.018 


0.082 


vt 


















1.000 


-0.653 


vt 




















1.000 



Table 4: Statistical correlations between selected parameters for Xp > 0.2. 
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Xp > u.o 


--y \ n /I 


Xp > u.o 


Xp > u.o 


Vd 


0.3866±0.0258 


0.19z4±0.01o0 


A AOOO 1 A AA/^A 

0.0888±0.0062 


A A /I AA 1 A AA01 

0.0400±0.0031 
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Vn 


U.oooi±U.Uzo ^ 
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0.0011±0.0047 
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T r-\- 


0.1227±0.0280 
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0.0376±0.0157 
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K± 
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n nni n_)_n nnAc, 

U .UU X U_1_U . UUtiO 


/ D 


0.0623±0.0069 


0.0190±0.0024 


0.0045±0.0009 


0.0007±0.0004 


0.0000±0.0008 


P 

Vd 


0.0362±0.0098 


0.0281±0.0100 


0.0056±0.0036 


0.0008±0.0023 


0.0000±0.0085 




0.0569±0.0139 


0.0183±0.0120 


0.0167±0.0058 


0.0050±0.0036 


0.0010±0.0055 




0.0328±0.0094 


0.0171±0.0059 


0.0051±0.0032 


0.0044±0.0021 


0.0014±0.0022 




0.0246±0.0092 


0.0076±0.0047 


0.0047±0.0032 


0.0006±0.0014 


0.0008±0.0033 


vl 


0.0232±0.0046 


0.0067±0.0015 


0.0016±0.0005 


0.0004±0.0002 


0.0002±0.0001 


Kg 


0.0461±0.0106 


0.0271±0.0144 


0.0193±0.0065 


0.0114±0.0034 


0.0041±0.0162 


Kg 
Tin 


0.0228±0.0146 


0.0126±0.0177 
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0.0000±0.0034 
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Kg 


0.1210±0.0100 


0.0743±0.0059 


0.0402±0.0040 


0.0192±0.0022 


0.0103±0.0020 




0.0586±0.0151 


0.0289±0.0081 


0.0099±0.00.55 


0.0013±0.0011 


0.0009±0.0010 


% 


0.0271±0.0025 


0.0088±0.0013 


0.0020±0.000.5 


0.0002±0.0003 


O.OOOOiO.OOOl 


^d^ 


0.0231±0.0032 


0.0074±0.0025 


0.0016±0.0022 


0.0017±0.0017 


0.0000±0.0003 


^u^ 


0.0211±0.0031 


0.0071±0.0024 


0.0016±0.0022 


0.0017±0.0017 


0.0000±0.0003 


^s^ 


0.0493±0.0062 


0.0240±0.0058 


0.0137±0.0056 


0.0077±0.0044 


0.0013±0.0043 




0.0295±0.0098 


0.0330±0.0088 


0.0162±0.0078 


0.0010±0.0036 


0.0022±0.0093 




0.0180±0.0016 


0.0045±0.0009 


0.0004±0.0009 


O.OOOOiO.OOOl 


O.OOOOiO.OOOl 




40.8 


40.6 


53.4 


39.4 


39.0 



Table 5: Results for various values of the minimum Xp requirement with statistical and 
systematic errors combined. In the last row, the of the solution is given. The number 
of parameters in the fit is 45, of which 32 are free, nine are fixed to their Monte Carlo 
values and four by hadronisation symmetries. 
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Source of Error 


vf 














7s(K±) 


Charged purity and eff. 


0.0254 


0.0114 


0.0226 


0.0176 


0.0106 


0.0010 


0.0004 


0.0615 


Kg purity and eff. 


0.0002 


0.0006 


0.0048 


0.0043 


0.0008 


0.0006 


0.0004 


0.0298 


A purity and eff. 


O.OOOf 


0.0001 


0.0005 


0.0005 


0.0007 


0.0029 


0.0017 


0.0011 


Charm tag purity and eff. 






0.0001 








0.0001 


0.0003 


Hemisphere correlations 


0.0012 


0.0049 


0.0050 


0.0033 


0.0024 


0.0024 


0.0057 


0.0293 


Other background 


0.0002 


0.0003 


0.0008 


0.0002 


0.0002 


0.0003 


0.0001 


0.0032 


Charm tag background 


0.0017 


0.0061 


0.0033 


0.0028 


0.0006 


0.0011 


0.0019 


0.0121 


Fixed quantities 


0.0002 


0.0008 


0.0006 


0.0003 


0.0001 


0.0001 


0.0003 


0.0026 


Hadronisation symmetries 


0.0033 


0.0002 


0.0039 


0.0016 


0.0005 


0.0004 


0.0005 


0.0188 


SRc 


0.0009 


0.0004 


0.0006 


0.0008 


0.0001 


0.0001 


0.0002 


0.0043 




0.0001 


0.0001 


0.0001 


0.0001 






0.0001 


0.0005 


Total Systematic Error 


0.0257 


0.0139 


0.0244 


0.0187 


0.0109 


0.0041 


0.0063 


0.0779 



Table 6: Systematic errors on the measurements of the 77*, corrected for detector efficiency 
and particle misassignment, for Xp > 0.2. Also shown in the last column are the systematic 
error contributions for 7s (K^). Absence of a number means that the error was less than 
5 X 10"^ 
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•^cut 






Xp > 0.2 


0.513±0.060±0.078 


0.381±0.069±0.034 


Xp > 0.3 


0.448±0.117±0.123 


0.365±0.127±0.152 


Xp > OA 


0.166±0.110±0.111 


0.180±0.111±0.087 


Xp > 0.5 


0.506±0.172±0.140 


0.593±0.140±0.160 


Xp > 0.2 


0.422 ± 0.049 ± 0.059 



Table 7: Results for different values of the minimum Xp cut for strange quark suppression 

as estimated by 7s(K^) = rj^'^/rif^ and 7s(Kg) = r]f^/r]^^ with statistical and 
systematic errors calculated taking into account correlations between the numerators and 
denominators of the ratios. In the last row a combined charged and neutral kaon result 
is given, corrected for decays. 









Xp > 0.2 


0.637±0.173±0.083 


0.468±0.069±0.030 


Xp > 0.3 


1.51 ±1.18 ±0.558 


0.307±0.112±0.046 


Xp > 0.4 


0.335±0.234±0.109 


0.118±0.141±0.076 


Xp > 0.5 


0.165±0.421±0.257 


0.221±0.264±0.173 



Table 8: Results for different values of the minimum Xp cut for uu diquark and strange 
diquark suppression as estimated by rj^/r]^ and rj^/r)^, respectively, with statistical and 
systematic errors calculated taking into account correlations between the numerators and 
denominators of the ratios. 



•^cut 


</< 




Xp > 0.2 


0.149±0.023±0.017 


0.184±0.030±0.007 


Xp > 0.3 


0.096±0.056±0.030 


0.162±0.041±0.006 


Xp > 0.4 


0.188±0.052±0.016 


0.170±0.071±0.014 


Xp > 0.5 


0.125±0.078±0.030 


0.200±0.117±0.013 



Table 9: Results for different values of the minimum Xp cut for ud diquark suppression as 
estimated by r]^/i]^ and with statistical and systematic errors calculated taking 

into account correlations between the numerators and denominators of the ratios. 
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Figure 1: Tagging probabilities as a function of the minimum Xp cut for charged pions 
and kaons. Data points are correlated for different values of the minimum Xp cut. The 
errors shown are statistical plus systematic. The lines show the JETSET predictions. 
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Figure 2: Tagging probabilities as a function of the minimum Xp cut for charged pions 
and kaons. Data points are correlated for different values of the minimum Xp cut. The 
errors shown are statistical plus systematic. The lines show the HERWIG predictions. 
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Figure 3: Tagging probabilities as a function of the minimum Xp cut for protons. Data 
points are correlated for different values of the minimum Xp cut. The errors shown are 
statistical plus systematic. The lines show the JETSET and HERWIG predictions. 
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Figure 4: Tagging probabilities as a function of the minimum Xp cut for A baryons. Data 
points are correlated for different values of the minimum Xp cut. The errors shown are 
statistical plus systematic. The lines show the JETSET and HERWIG predictions. 



30 



^ 25 



^ 20 
a 

> 

W 15 



10 
5 





□_ 

1.08 1.1 1.12 1.14 

Invariant mass pit (GeV) 



1.1 1.12 1.14 
Invariant mass p7U (GeV) 



> 


10 


Me 


9 


go 


8 






a 


7 






> 


6 




5 




4 




3 




2 




1 








T 



n — I — I — I — I — I — r 



T 



AD , AD 
Xp>0.5 



5 

- 1/1 



_L 



35 
30 



1.08 1.1 1.12 1.14 

Invariant mass pTi (GeV) 



a 25 
W 20 



15 
10 
5 




T 1 1 1 1 1 T 



AK^, AK 
Xp>0.5 



T 



T 




1.08 1.1 1.12 1.14 

Invariant mass pTi (GeV) 



Figure 5: Invariant mass distributions for the A signals in data (points with error bars) 
and the JETSET Monte Carlo (histogram). The Monte Carlo is normalised to the same 
number of events with a D*~ (left) or a K"*" (right) in the opposite hemisphere, for two 
different Xp ranges. 
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Figure 6: Determination of 7s = s/u using the estimator rj^ /rjf (a) and Vd^/v^^ 
the (b) in the OPAL data (sohd points). Data points are correlated for different values 
of the minimum Xp cut and the errors shown are statistical plus systematic. The solid 

lines represent the true "rj^^ /rjf^ and f]^^ /r^^ in the JETSET Monte Carlo and the 
dashed lines the HERWIG predictions. The dotted lines represent the input value of 
7s = PARJ(2) = 0.31 in JETSET. 
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Figure 7: Vd/Vu (^) Vd/Vs 0^) OPAL data (solid points). Data points are 

correlated for different values of the minimum Xp cut and the errors shown are statistical 
plus systematic. The solid lines represent the true 't]'^/'!]^ and Tj^/rj^ in the JETSET 
Monte Carlo and the dashed lines the HERWIG predictions. The dotted lines represent 
the input values of uu/ud = 3 • PARJ(4) = 0.075 (a) and us/ud = 7^ • PARJ(3) = 0.1395 
(b) in JETSET. 
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Figure 8: Vu/Vu (^) ^^'^ vt/vf' (b) in the OPAL data (solid points). Data points are 
correlated for different values of the minimum Xp cut and the errors shown are statistical 
plus systematic. The sohd lines represent the true r]^/r]^ and ■q^/f]^ in the JETSET 
Monte Carlo and the dashed lines the HERWIG predictions. The dotted hues show the 
input values for qq/q = PARJ(l) = 0.085 in JETSET. 
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